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ABSTRACT 

In a weakly ionized plasma, the evolution of the magnetic field is described 
00 , by a "generalized Ohm's law" that includes the Hall effect and the ambipo- 

lar diffusion terms. These terms introduce additional spatial and time scales 

m 

' which play a decisive role in the cascading and the dissipation mechanisms in 

| magnetohydrodynamic turbulence. We determine the Kolmogorov dissipation 

scales for the viscous, the resistive and the ambipolar dissipation mechanisms. 

•rH . 

■ The plasma, depending on its properties and the energy injection rate, may 

preferentially select one of the these dissipation scales, thus determining the 
shortest spatial scale of the supposedly self-similar spectral distribution of the 
magnetic field. The results are illustrated taking the partially ionized part of 
the solar atmosphere as an example. Thus the shortest spatial scale of the 
supposedly self-similar spectral distribution of the solar magnetic field is de- 
termined by any of the four dissipation scales given by the viscosity, the Spizer 
resistivity ( electron-ion collisions), the resistivity due to electron- neutral col- 
lisions and the ambipolar diffusivity. It is found that the ambipolar diffusion 
dominates for resonably large energy injection rate. The robustness of the 
magnetic helicity in the partially ionized solar atmosphere would facilitate 
the formation of self-organized vortical structures. 

Key words: Partially ionized plasma, Kolmogorov dissipation, Hall effect, 
Ambipolar diffusion. 
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1 INTRODUCTION 

It is widely recognized that some ideal integrals of motion (i.e., constants of motion in the 
dissipation-less limit) play an essential role in the complex nonlinear dynamics leading to 
remarkable self-organizing structures. The "helicity" is one such quantity that imposes a 
"topological" constraint on a field Moffat (1978). Constancy of the helicity is not attributed 
to any "symmetry" of the system, but is due to a "defect" (or a singularity) of the Poisson- 
bracket operator, and, thus, is a robust constraint throughout the evolution. 

A dissipative mechanism can change the helicity to remove the topological constraint. 
Self-organization may be understood as a subtle balance between the conservation (restric- 
tion) and the dissipation (relaxation). 

One of the most successful models of self-organization is due to Taylor (1974), who 
invoked the constancy of the magnetic helicity H = J A- B d 3 x/2 (the integral is taken over 
the total volume). Minimizing the magnetic energy E — J \B\ 2 d 3 x/2 (he omitted the kinetic 
and the thermal energies) for a fixed H yields an Euler-Lagrange equation VxB = XB (A is 
the Lagrange multiplier for restricting H) , whose solution is the "Beltrami field" representing 
a twisted force-free magnetic field. The cross helicity H c — JV-B d 3 x/2 is also an invariant 
of the ideal magnetohydrodynamic system with flows where V is the fluid velocity. 

This "variational principle" is based on the assumption that the energy E is preferentially 
(selectively) dissipated while the helicity H is approximately conserved (and that no other 
conserved quantity puts an obstacle for the minimization of E). A possible justification 
for this assumption is given by the "energy cascade" ansatz — the energy density of the 
fluctuating field cascades toward small spatial scales, and the variations (spatial derivatives) 
of the field are enhanced. Then, the dissipation of the energy, which includes higher-order 
spatial derivatives in comparison with the helicity, proceeds much faster than that of the 
helicity. This "scale change" of fluctuations may be represented by the cascade of the energy 
density in the Fourier (wave number) space. 

It has been established that E and H are the integrals of motion of an ideal magnetohy- 
drodynamic system, the energy cascade yields the selective dissipation of the energy (with 
respect to the helicity) leading to the Taylor relaxed state. This paper extends the scope 
of these considerations to "weakly ionized plasmas" . There are many astrophysical systems 
with a rather low degree of ionization dominated by the charged particle-neutral collisions 
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and the neutral particle dynamics. A major part of the solar photosphere (Leake and Arber 
2006; Krishan and Varghese 2008;), the protoplanetary disks (Krishan and Yoshida 2006) and 
the molecular clouds (Brandenburg and Zweibel 1994) are some of the examples of weakly 
ionized astrophysical plasmas. In addition these systems are believed to be turbulent. The 
evolution of the magnetic fields in such a plasma would be affected by the multifluid inter- 
actions in general and the ambipolar diffusion in particular (Zweibel 1988). It is important 
to know which dissipation scale would be the most effective in a given situation. 

Here, We (1) set up the energy and the helicity evolution equations in a weakly ion- 
ized plasma and (2) determine the Kolmogorov dissipation scales for the different energy 
dissipation mechanisms. The weakly ionized part of the solar atmosphere is presented as an 
example to illustrate the predominance of one or the other dissipation scale. We show that 
the ambipolar effect, which is nonlinearly enhanced when the energy injection rate is large, 
may dominate the energy dissipation, while it conserves the helicity. 



2 ENERGY AND HELICITIES IN WEAKLY IONIZED PLASMAS 

The dynamics of a weakly ionized plasma can be described with the following equations ( 
Krishan and Yoshida 2006, Krishan and Varghese 2008, Krishan and Gangadhara, 2008 ): 
the momentum balance of the neutral fluid is given by: 

^ + (v ■ v) v = -vh + + /iW, (i) 

at cp n 
and the magnetic field B evolves as: 



dB „ 

sr = Vx 

and 



J JxB 

V + 

en e cUi n pi 



xB\+r]V 2 B. (2) 



47T 

VxB = ™J (3) 

where V is the velocity of the neutral fluid and h is the total enthalpy. The Lorentz force in 
the neutral dynamics appears due to the ion- neutral coupling in that the Lorentz force on 
the ions (erii(E + Vi * B )) is balanced by the ion-neutral collisional force (— PiVin{Vi — V n )) 
where u in is the ion- neutral collision frequency. The ion inertial force has been neglected 
in comparison with the ion-neutral frictional force. The viscosity \x of the neutral fluid 
arises due to neutral- neutral collisions with frequency v nn . The electron inertial force is also 
neglected. We have further assumed a constant density incompressible system. We may write 
the electric diffusivity r\ as: 
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2 

C Ul e V en o i A \ 

v= tol&r = s ° Vm > (4) 

and the kinematic viscosity p as: 

v- = — = K v ™» ( 5 ) 



where 5 e := c/uj pe = c/ \J 4irn e e 2 j 'm e is the electron inertial length, and A n := v n) th/v nn is the 
neutral mean free path: u en , v nn and v e i are the electron-neutral, the neutral-neutral, and 
the electron- ion collision frequencies; v e i << v en in the case of the weakly ionized plasma. 

Assuming constant densities and incompressibility, we normalize the variables in the 
'Alfvenic units" . We chose a representative neutral flow velocity Vq and a magnetic field Bq 
related as: 

„ t/2 r>2 

2 8tt' 1 ' 

Solving (jHJ) for Vq gives a virtual Alfvenic velocity Vq = Bq/ (Ait p n ) 1 / 2 . If p n were pi, this Vq 
is the well-known Alfven velocity. In a weakly ionized plasma, 

a := ^ > 1 (7) 
Pi 

For example, in the solar atmosphere, a can be as large as 10 6 , the protoplanetary disks 
have a of the order of 10 8 . 

Let Lq be a characteristic length in the system. We normalize x by L , t by to := Lq/Vq, 
B by S , V by Vb and energy densities by p n V 2 = Bq/(8tt). Using these variables, Equations 
(1) and (2) can be written in the (pQ) and ([2]) dimension-less form: 
dV 

— + (V- V)V = -Vh + JxB + e M V 2 V, (8) 

ap 

- = Vx{[V- fH J + fi Jx5]xB}| £r ,V 2 B, (9) 

where the scaling parameters are defined as: 

c/uj p i 5i 

e H ■ = a—^- = a—, (10) 
ga ■= e H — , (11) 



e 



to. I 

V 'I 7-2 ' T y 



^ = *tV. (12) 



-^o -^o ^0 

Here, 5« is the ion inertial length. The scaling parameter e# multiplying the Hall term is 
enhanced by the factor a = p n j pi in comparison with the standard (fully ionized) Hall term. 
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With an appropriate homogeneous boundary conditions, we have the energy equation 
(E = J(V 2 + B 2 )dx/2), 

^L = - eA J jIb* d 3 x -e v Jj 2 d 3 x - e M J | V x V\ 2 d 3 x, (14) 
and, for completeness, the helicity equation (H = J A ■ Bdx/2) 

^L = - e [j. B d 3 x. (15) 

In the ideal limit (ea = 0, e v = and e M = 0), E and H are conserved. The constancy 
of H is destroyed only by a finite resistivity e v . The cross helicity H c of the ideal MHD 
transforms to the ion canonical helicity Hq = J(A + e^F) ■ V x (A + EnV)dx/2. The cross 
helicity H c is conserved when e# = 0, €a = and the ion helicity H G is conserved for 6^ = 0. 
For 6a 7^ both if c and ifc are not conserved. From ( fl4|) . we find that the energy dissipation 
is contributed by (i) ambipolar diffusion (scaled by e^), (ii) resistivity (friction of electrons 
with neutrals; scaled by e^), and (iii) neutral viscosity (scaled by e M ).The advective terms 
(V ■ V)V, V x (V x B) and V x (J x B) ( the Hall term ) are responsible for the energy 
cascade mechanism. Since the velocity and the magnetic field are coupled through equations 
(8) and (9), all the advective and the dissipative processes operate on both the velocity 
and the magnetic fields. We shall define the Reynolds number in a rather broad sense as 
the ratio of the advective term and the dissipation term. We will see that the Reynolds 
numbers resulting from different combinations of the advective process and the dissipation 
process set up a rather complex scale hierarchy of the Kolmogorov microscale, in the wave 
number (k) space, which changes depending on the strength of the energy injection rate. 
In what follows we will select one dissipation mechanism and combine it with the three 
advective ( cascading) mechanisms to define the three Reynolds numbers and determine the 
corresponding Kolmogorov microscales. 



3 SCALE HIERARCHY 

The wave-number (k) space is divided into (i) the energy injection (large scale) range, (ii) 
the "inertial range" , and (iii) the dissipation range. The inertial range is dominated by the 
convective [(V ■ V)V], inductive [V x (V x B)] and the Hall [-V x (e H J x B)] effects, 
which create a sub-hierarchy in the inertial range. The higher-fc-end of the inertial range 
is the "Kolmogorov scale" that is determined by one of the three (ambipolar, viscous and 
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resistive) dissipation mechanisms. The aim of this section is to estimate the Kolmogorov 
scale by identifying the responsible (i.e., the dominant) mechanism of energy dissipation. 

Here we invoke Kolmogorov's ansatz of "local interactions" in the /c-space, which may 
be formulated as follows. Let K symbolize a "range" of wave-number space, i.e., K stands 
for {k; K ^ \k\ < K + A}, where k is the wave vector and A is a certain positive constant. 
Denoting by u(k) the Fourier transform of a field u{x), we define u K = J2keK u(k)e' lk x , which 
means the component of u(x) in the hierarchy of wave-numbers ranging in K. Suppose a 
term X is included in an evolution equation of a physical quantity u. When we observe the 
hierarchy K, Xk contributes the temporal variation (du/dt)x- If X is a linear term including 
a field v (and, possibly, the differential operator V), then we may estimate X K only by vk- 
But, if X includes, for example, v ■ w, then all v(ki) and w(k 2 ) (and the corresponding V 
translating into iki and ik 2 ) may contribute to X K , if ki ± k 2 £ K. Now, the "locality 
ansatz" claims that only vk and wk (and, thus, V of order K) dominates Xk- We assume 
that this ansatz holds for all (even higher-order) nonlinear terms in our system. We also 
assume that no geometric anisotropy diminishes the magnitudes of vector products such as 
J x B. Hence, we estimate, for example, | J x B\ K « | J\ K \B\ K « K\B\ 2 K (using J = V x B 
in the normalized unit). In what follows, we denote \B\ K by B K . 

3.1 Kolmogorov scales defined by viscosity dissipation 

Using a rather general definition of the Reynolds number as the ratio of the term responsible 
for the cascade process to that responsible for the dissipation process, we consider, first, the 
advective term (V ■ V)V and the dissipation due to the viscosity. This defines the very 
familiar, the (standard) Kolmogorov scale, when the other mechanisms (ambipolar diffusion 
and resistivity) are still negligible. 
The Reynolds number is defined by 



|(V-V)V 


\K 






K 



(16) 



Note that this Reynolds number is evaluated for each scale hierarchy as a function of K. 
Invoking the previous assumptions, we estimate 

R,(K) « \ (17) 

The "Kolmogorov microscale" is usually defined as the spatial scale / at which the Reynolds 
number becomes unity ( e.g. Holmes, Lumley and Berkoz 1996, Frisch 1995 ). We prefer 
to work in the wavevector space such that K = l^ 1 and define the inverse Kolmogorov 
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microscale K = where R^(K) becomes of the order of unity ( e.g. Pope 2000). We 
obtain 

V K 

K, = (18) 

To estimate (eliminate) Vr , we invoke the energy cascade rate ( total energy E being an 
invariant of the system) £ ( normalized by So = V 3 j 'Lq ) that is assumed to be scale invariant 
and equal to the energy dissipation rate as well as the energy injection rate. From Eq.(13) 

£ = ^K V l = (is) 

Solving ffl9l for Vjf , we obtain 

V K(l = £ l ^ A . (20) 
Substituting (I2D1 into ({TBI yields the well-known relation 

k, = ev%-w (21) 

If the energy cascade is dominated by the advective term V x (V x B) and the dissipation 
by the viscosity; the Reynolds number R^ M (K) should better represent the relation between 

the energy cascade and the dissipation where 
\Vx(VxB)\ K ^ B K V K 
e^\V 2 V\ K e^K e^KCy/B 

where we have introduced a coefficient Cy/B '■= Vk/Bk that is, in general, a function of K. 
The corresponding Kolmogorov scale becomes 

Bk^m VlCjjM cl/i^ -3/4^-1 / O o\ 

A m m - — ^ - t> e M o v/B . [26) 

In the ideal MHD regime, however, we may assume Cy/B ~ 1 throughout the MHD iner- 

tial range, because the energy transfer is an equal collaboration of the induction and the 

convection terms. Hence, K^m ~ K^. 

The situation changes, when the Hall term dominates the energy cascade. This is the 

case if the ion skin depth multiplied by the density ratio a, en (normalized by the system 

size) is larger than the Kolmogorov length scale, i.e., e#i^ > 1. One can define the Hall 

Reynolds number by taking the advective term e# J x B ( the Hall term) and the dissipative 

term to be due to the viscosity. The corresponding Reynolds number is defined as 

p (v\ kffV x (J x B)\k _ e H B 2 K 6hVk , . 

H^ H [K):= Tv~ = e ^ ' 

CmI V V \K tfiVK e ^V/B 

scales the ratio of the energy cascade rate and the dissipation. Assuming Cy/B = (^hK) p 
(for K > 1/en) with a certain exponent p (the Hall-MHD turbulence theory, Krishan and 
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Mahajan 2005, predicts, for R-mode turbulence p = 1, while for L-mode turbulence p = — 1), 
we obtain a "Hall Kolmogorov scale" (denoting q := l/(2p)) 

K,n = elT 1 = (^'V 3 ' 4 )* = (e^p (K.r • (25) 

Since K^ H /K^ = l/^enK^) 1-9 (and we are assuming e^if^ > 1), K^ H interpolates ejj 1 and 
(i.e., e^ 1 < K^ H < K^) as long as q < 1. 

As mentioned above, the condition for the appearance of the Hall MHD regime is 
e H K, = ene^E 1 " > 1. (26) 

3.2 Kolmogorov scales defined by resistivity dissipation 

Next, we examine the case where the resistive dissipation dominates over the viscous and 
the ambipolar diffusivities. Both the viscous and the resistive dissipations have a common 
mathematical structure viz. they are linear terms multiplied by V 2 . 

The resistive and the viscous dissipation mechanisms may be compared as follows: from 
(114p . the ratio of the resistive and viscous dissipation terms is 

|e„V x B\ 2 K _ t n K 2 B 2 K _ fv_(j-2 ^27) 



|e M VxV|! H KW 2 e^ v ' B - 

Again taking (V ■ V) V) as the advective term, in the MHD regime where we may assume 
Cv/b ~ 1, the resistive Kolmogorov scale is given by just replacing e M by e„ in (12T1) . i.e., 

K v = £ 1 '%~V\ (28) 

where the energy dissipation rate is £ = e^K^B"^ . Evaluating 8 for the same energy injection 
rate, we may write K v = (e„/e M )~ 3 / 4 iT M . By (0J and ([5]), we estimate 

ej, = |^en (2Q) 

For the MHD regime where we may assume Cy/B ~ 1 the advective term V x (V x B) 
along with the resistive dissipation furnishes the same scale K v . 

In the Hall MHD regime {ehK v > 1), we replace the previous Hall Reynolds number 
iBb by 



R M K) : = \Sl " M * B ^ „ Z**. (30) 
\e v v t»\ K e v 

Hence, at the corresponding Kolmogorov scale K v h, we estimate Bk H = e„/e#. From the 
energy dissipation rate £ = e n K 2 H B\ H , we obtain 

K vH = £ 1/2 e~ 3/2 e H = K 2 e H . (31) 
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Hence, unlike the relation of the viscosity Kolmogorov scales ( 1251) . K v h > K^. 



3.3 Kolmogorov scales defined by ambipolar diffusion 

Finally, we estimate the Kolmogorov scales assuming that the dissipation is dominated by 
the ambipolar term. 

We begin with the MHD regime where we may assume Cy/s ~ 1. As shown earlier in this 
case the advective terms (V ■ V) V) and V x (V x B) furnish equal Kolmogorov microscales. 
We define an "ambipolar Reynolds number" by 
RJK) .- Wx(VxB)\ K Vk 

Ka[K) - \e A V x [(J x B) x B]\ K ~ e A KB%' {6l) 
The ambipolar Kolmogorov scale K A is characterized by R A (K A ) « 1: 

= = fW- < 33 > 

The ambipolar dissipation rate is now ( Eq. 13): 

S = e A K\B± KA = e- A x C vjB B\ A . (34) 
Plugging B Ka = £ 1 / 2 e A 2 C v j B into fl33l . and assuming Cy/B ~ 1, we obtain 
K A = e A 3/2 £^. (35) 
Let us compare (|35|) with (I2T]) : 

*j± = efe-*'*£-*l\ (36) 

From this relation, we see that the ambipolar diffusion dominates over the viscous dissipation 
(i.e., K A j Kfj, <C 1), when the energy dissipation rate (=energy injection rate) S is sufficiently 
large. This is because the ambipolar dissipation is a nonlinear term that is enhanced when 
the fluctuation level is high. 

If the Hall term dominates the energy cascade i.e. taking e#V x ( J x B as the advective 
term (euK A > 1), we replace f[3"2j) by 

Rah{K) - \e A V x [(J x B) x B]\ K ~ 7^- K (3?) 
Using fl34|) . we estimate -Ba'ah = (^/ e A) 1 ^ 4 K A H^ which, together with R AH (K AH ) = 1 yields 

K AH = eTe H 2 £ 1 ' 2 = {e H K A )~ 2 K A . (38) 

We, thus, see K AH < K A . 

Comparing with (125]) . we get: 
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dissipation mechanism MHD regime (ejjif < 1) H MHD regime (e[fK > 1) 

viscosity K„ = £ l / 4 e~ 3/4 = (e^ 1 ) 1 -^* (< K„) 

resistivity K v = £ 1 ' i e~' i/i K vH = e H K 2 (> K v ) 

ambipolar K A = £' 1 / 2 e~ A 3/2 K AH =e H 2 R- 1 « K A ) 

Table 1 . Summary of Kolmogorov scales 



Kah = e 3/2 e -(g+l) e 3g/4£;(2-g)/4 = ( ^ei\ 3/2 £ (l/2)-? e 3f/4£(2-g)/4_ ^ 



3.4 Comparison of different Kolmogorov scales 

In Table UJ we summarize different estimates of Kolmogorov scales. When the plasma pa- 
rameters and the energy injection rate £ are specified, we can estimate the appropriate 
Kolmogorov scale determined by the balance between the relevant energy cascade and the 
dissipation mechanisms. 



4 SCALE HIERARCHY ON THE SOLAR ATMOSPHERE 

The solar magnetic flux, believed to be generated in the convection zone, has to pass through 
the partially ionized solar photosphere before it can appear high up in the solar corona. This 
realization is rather recent and is now receiving a lot of attention. Arber, Haynes and Leake 
(2007) has emphasized the profound effects on the temperature and the current structure 
of the overlying chromosphere and the corona that the inclusion of the neutral medium can 
produce. It is important to know which dissipation scale would be the most effective for the 
conditions typical of the solar atmosphere. We estimate the various Kolmogorov scales in 
the partially ionized part of the solar atmosphere. 

The various collision frequencies are determined from Khodochenko et al.(2004): 

(8K T \ 
— > (40) 
Trmy / 

where (i,j) stands for the species of particles, the cross-section S e „ w 10 _15 cm 2 , Sj n w 
E nn ~ 5 x 10~ 15 cm 2 , rriij = mivrij (m^ + m^) -1 , n n is the neutral particle density and T is 
temperature in degree Kelvin. The typical values of the physical parameters on the solar 
atmosphere are given in Table II. 

The variation of the ionization fraction with height on the solar atmosphere is shown 
in figure (1). 

The ratios of the different Kolmogorov scales can be expressed as: 



Kolmogorov Dissipation scales in Weakly Ionized Plasmas 11 



h 


T 


Pi 




Pn 




B 


Be 




0. 


6520. 


1 no * 1 n~ 


10 


1 Qn v 1 n~ 

J- . 3U A IU 


07 


1 9nn no 


9.4x 10 


-04 


50. 


5790. 


1 90vin~ 

j_ . z*yj A iu 


ii 


1 ^QvlfP 
± . jy a ±u 


-07 


1125.77 


1.2x 10 


-04 


±zo. 


OZ 1 u. 


1 isy1 n - 

l.loX 1U 


12 


l.UU X 1U 


-07 


you. id 


l.O X 1U 


-05 


175. 


5060. 


3.39x10" 


13 


7.04x10- 


08 


880.33 


4.6x10 


-06 


250. 


4880. 


9.37x10" 


14 


3.89x10- 


-08 


737.21 


1.7x10 


-06 


400. 


4560. 


1.12x10" 


14 


1.09x10- 


-08 


503.71 


4.2x10 


-07 


490. 


4410. 


4.37x10- 


IS 


4.84x10- 


09 


394.42 


2.6x10 


-07 


560. 


4430. 


4.72x10- 


15 


2.47x10- 


-09 


322.27 


4.2x10 


-07 


650. 


4750. 


2.29x10- 


14 


1.00x10- 


-09 


246.31 


3.7x10 


-06 


775. 


5280. 


1.08x10" 


13 


3.79x10- 


10 


183.67 


3.5x10 


-05 


855. 


5650. 


1.75x10" 


13 


1.66x10- 


10 


143.40 


1.0x10 


-04 


980. 


5900. 


1.78x10" 


13 


6.57x10- 


-11 


108.65 


1.8x10 


-04 


1065. 


6040. 


1.67x10" 


13 


3.60x10- 


-11 


90.88 


2.5x10' 


-04 




i(r 4 



lo-' 




600 
height 



1200 



Figure 1. Variation of the Ionization fraction with height h on the solar atmosphere. 



= 3 x 10 18 S~ 3 / 4 pTpn 3/4 T 9 / 8 B- 3 / 2 , 



Ka 

= io-W /4 p; 3/2 



(41) 
(42) 



We present plots of the ratios j^- and -^A vs height for different values of the energy- 
injection rate £ vs height in Figs. (2), (3) and (4) respectively. 

We observe from Fig. (2) that the ambipolar dissipation predominates over the viscous 
dissipation i.e. Ka < K^ on major part of the solar photosphere and the chromosphere. 
Choosing the parameters at a height of 980 Km above the photosphere as the normalizing 
parameters i.e. p n = 6.5 x 1CT 11 g cm -3 , B = 108 G, Lq ~ 1000 Km, we find Vq ~ 3.5 Km 
s _1 and £, = 4 x 10 8 cm 2 s~ 3 . The condition for the predominance of the ambipolar diffusion 
over the viscous dissipation becomes £ > 1.2 x 10 -2 or, recovering dimensions, 5 x 10 6 cm 2 s~ 3 . 
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Figure 2. Variation of the ratio with height h for different values of the energy injection rate £ 
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Figure 3. Variation of the ratio with height h 
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One can estimate the typical injection rate of the turbulent convective energy on the sun to 
be V^/Lq ks 10 8 cm 2 s" 3 by taking a typical convective velocity V c w 2Kms~ 1 . The condition 
for the predominance of the ambipolar diffusion is, thus, easily satisfied. The inclusion of 
the Hall effect modifies Ka to Kah and as seen from table I Kah < Ka, reconfirming the 
predominance of the ambipolar effect as the dissipation mechanism. Figure (3) demonstrates 
the predominance of the resistive dissipation (due to electron-neutral collisions) over the 
viscous dissipation. 
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Figure 4. Variation of the ratio ^-A with height h for different values of the energy injection rate E 

That the ambipolar dissipation predominates over the resistive dissipation i.e. Ka < K v 
for the region extending from a height of 200 Km to 1000 Km for reasonable values of 
the energy injection rate £ can be seen in figure (4). The Spitzer resistivity in a weakly 
ionized plasma, by definition, is much smaller than the resistivity due to the electron-neutral 
collisions. We conclude that the ambipolar dissipation would dominate over other dissipation 
mechanisms on most of the partially ionized region of the solar atmosphere. 

5 CONCLUSION 

The dynamics of weakly ionized plasmas is governed by the neutral fluid being subjected to 
the Lorentz force and the magnetic induction to the Hall and and the ambipolar effects. There 
are now three mechanisms by which the energy could dissipate, each mechanism having its 
own characteristic dissipation scale. We have determined these scales as presented in Table I. 
The example of the solar atmosphere shows that the ambipolar dissipation would determine 
the short scale end of the supposedly self-similar distribution of the photospheric magnetic 
field thereby identifying the region of dissipation leading to heating. We observe that H 
is a more robust quantity in comparison with E, if either the ambipolar or the viscous 
dissipation is larger than the resistive dissipation. In the standard argument of selective 
dissipation, moreover, even the resistivity dissipates E faster than it does H because of 
the energy cascade toward small scales; dE/dt includes higher-order spatial derivatives in 
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comparison with dH/dt, so that dE/dt assumes a larger value for small-scale fields. This 
could lead to the formation of organized plasma structures, vortical structures with twisted 
magnetic fields, so ubiquitous on the sun. 
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